Entanglement Generation via a Completely Mixed Nuclear Spin Bath 



H. Christ, J. I. Cirac, and G. Giedke 
Max-Planck-Institut fur Quantenoptik, Hans-Kopfermann-Str. 1, D-85748 Garching, Germany 

(Dated: September 13, 2008) 

We show that qubits coupled sequentially to a mesoscopic static completely mixed spin bath 
via the Heisenberg interaction can become highly entangled. Straightforward protocols for the 
generation of multipartite entangled (Greenberger-Horne-Zeilinger-)states are presented. We show 
the feasibility of an experimental realization in a quantum dot by the hyperfine interaction of an 
electron with the nuclear spins. 
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I. INTRODUCTION 

The quest to realize quantum information process- 
ing (QIP) has motivated an impressive race to imple- 
ment high precision preparation and manipulation of iso- 
lated two-level quantum systems (qubits) in a wide va- 
riety of physical settings^. A hallmark achievement for 
each such approach is the generation of quantum entan- 
glement through controlled interaction between two or 
more qubits. Since switchable direct interactions between 
qubits often entail additional decoherence mechanisms, 
many QIP proposals rely on interactions mediated by an 
additional quantum system. As a rule this mediator (just 
as the qubits themselves) needs to be prepared in a pure 
state to achieve high-fidelity quantum operations and it 
may look futile to use a high-entropy mesoscopic spin 
bath for this task. In contrast to these expectations, we 
show here that high-fidelity entanglement generation can 
be realized even if the qubits can only interact with an 
arbitrarily mixed spin bath, provided that this interac- 
tion can be switched on and off, single-qubit unitaries 
are available, and the bath has slow internal dynamics. 
This is motivated by and will be illustrated through the 
example of electron-spin qubits in quantum dots (QDs) 2 , 
where the ensemble of lattice nuclear spins represents a 
strongly coupled but slowly evolving spin bath. 

Nuclear spins in quantum dots have received much 
theoretica l 3 * 4 * 5 * 6 * ^ * and experimental^^ attention in the 
QIP context as the main source of electron-spin deco- 
herence through the strong hyperfine coupling. It has 
also been noted that their slow internal dynamics and 
long (expected) decoherence time^ make the ensemble 
of nuclear spins useful as a quantum memory® or for 
quantum computation^. These applications, however, 
require careful and yet unachieved preparation of the nu- 
clear system. What we show here is, that the unpre- 
pared, highly or even maximally mixed (nuclear) sys- 
tem is able to mediate coherent interaction between elec- 
trons and thereby allows the generation of highly entan- 
gled states of many (electron spin) qubits without any 
electron-electron interaction. 

We consider a QD in the single-electron regime^ and 
assume the availability of single-electron state prepara- 
tion and measurement as well as the controlled shuttling 



of prepared electrons into and out of the QD, all of which 
have been demonstrated experimentally^. Additionally 
required is control of the detuning (e.g., by a magnetic or 
electric field), which switches the hyperfine (HF) interac- 
tion between resonant and off-resonant regimes. We first 
show how sequential interaction of three electrons with 
the nuclear bath can generate a maximally entangled pair 
of electron spins. More generally, the class of states that 
can be generated via the spin bath is characterized in 
terms of matrix product states. Finally we show that 
imperfect electron spin operations, inhomogeneous cou- 
plings between electron and nuclei and modifications to 
the ideal static spin bath still allow for the scheme to be 
realized. In situations where spin-orbit coupling is large, 
our scheme can be an interesting alternative to the stan- 
dard exchange based setups, because it doe s not involve 
occupation of any higher orbital levels^^*. 

II. ENTANGLEMENT GENERATION 

We consider each electron coupled via the uniform 
Heisenberg interaction to the bath of N nuclear spins 
and to an external magnetic field B z (h = 1) 

H= 4t + S+I ^ + P* 8 * + 9>bB,S*. (1) 

S is the spin operator for the electron and I M — J2i 
are the three components of the collective nuclear spin 
operators (n = ±,z and [7+, I z ] = [I + ,I~] = 2F). 

g* is the electron ^-factor and [Ib the Bohr magneton. 
We consider spin-1/2 nuclei, neglect bath dynamics, the 
bath spins' Zeeman energies, and inhomogeneities in the 
Heisenberg couplings for now, and discuss the validity of 
these approximations towards the end of this article. 

We use the Dicke basis {| I, m, /?)}, where 1(1 + 1) is 
the eigenvalue of the collective angular momentum oper- 
ator I 2 , the eigenvalue of I z isgiven by m, and (3 is the 
permutation quantum number^S. The initial state of the 
spin bath in the following is the identity 

Pbath - ^ E I 7 ' m ' ^X 7 ' TO > P\ = l2«/2 W (2) 

In the following we omit j3, which does not enter in the 
dynamics. This situation of a completely unknown bath 
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state is, e.g., a suitable description for GaAs QDs even 
at temperatures as low as 100 mKP^l. In the following, 
time will be given in units of N/A. Even though the 
idea we present is applicable to any (quasi-)static bath, 
we perform all estimations for GaAs, i.e. in particular 
A- 1 w 40 ps. 
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FIG. 1: (Color online) Left: Sketch of the protocol, (a) The 
z-polarized "control electron" interacts resonantly with the 
nuclear spin bath, (b) A sequence of a;-polarized electrons 
interacts off-resonantly with the bath, (c) The control elec- 
tron interacts resonantly again and is then measured in the 
z-basis. Right: Time dependence of overlap F with Bell state 
\4>-) = (| + +) - | - -))/V2 (solid blue line) for N = 10 3 . 
Dashed, red line shows the probability P for a [ measurement. 

The first electron spin (which we also refer to as an- 
cilla electron) is prepared in the state ||) and interacts 
resonantly for a time t\ with the nuclear spin bath 

U\I,m,V = c Im {ti)\I,m,1) + s Im (ti)\I,m+l,l), (3) 

with U = er lHtl and 

.'(l + 2J)iA ,l + 2m . /(l + 2i>i 
ci m (h) = cos I ; I - i- — sin 
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Then the next electron spin, initial state |+) = 
l/V2(|t) + |l)), interacts for 

a time ^2 off-resonantly (e.g., 
in the presence of a large B z ) with the spin bath. For 
g* /j,bB z 3> A/vN, the flip-flop part of the Hamiltonian 
can be approximately neglected^, yielding 
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where B = g* fisB z N/A and the index / has been omit- 
ted for brevity. Remarkably, by choosing the interaction 
time t% = it, the state of the electron is transformed to 
{-i) m |(-) m ), i.e. for even m = 2k to (-l) fc |+) and for 
odd m = 2k + 1 to — i{— l) fe | — ). For convenience we as- 
sume Btill = 2tt£, feN, which is adjusted by the "free" 
parameter of the large field. 

With the third electron, also in |+) initially and the 
same interaction, the state becomes 

±c Im (t 1 )\I,m,V\±±)Tsi m (ti)\I,m+ 1,|)|++), (4) 



with upper (lower) signs referring to even (odd) m. 

In the final step, the ancilla electron interacts reso- 
nantly with the nuclei again [cf. Eq. giving 

±c Jm (*i)| ±±)[c Jm (ti)|m,t) +si m {ti)\m+l,l)] (5) 
+s/m(*i)| TT)[c* Im (ti)\m+ 1,J.) + si m (ti)\m, T)], 

for even/odd m and is eventually measured projectively 
in the z-basis. If the measurement outcome is J, it is clear 
from Eq. ([5| that in each subspace the second and third 
electrons are in the maximally entangled state 
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where the phase <f) m — 2 arg(cj m ) depends on the quan- 
tum numbers / and m, leading to a washing out of the 
entanglement when the average over the different sub- 
spaces is taken. However, for short times ((2/+ l)ti <C 1 
for typical values of / ~ yN) this phase tends to zero 
and near ideal entanglement is created, albeit at the price 
of a lower success probability, see Fig. [T] 



III. MULTIPARTITE ENTANGLEMENT 

The presented scheme generalizes in a straightforward 
manner to multipartite entanglement creation. Following 
the same protocol using n electrons with arbitrary initial 
states • ■ ■ , \ip n ), the final state becomes 

|*„) - 1/V2 (1 + [(-l) m+ ha z f n ) |^i, • . ■ , Vn), (7) 

where the matrices are given in the standard z-basis and 
we assumed the short time limit t\ — + for clarity. If 
1^.) = |+) for all k, this is a n-partite Greenberger- 
Horne-Zeilinger (GHZ) state. The m-dependent relative 
phase in the above equation restricts to generation of 
GHZ-states with even particle number. 

When multiple resonant interactions with the ancilla 
and varying interaction times are allowed, a larger class 
of states becomes accessible. To see which states can the 
in principle be prepared, we exploit the similarity of our 
setup to the sequential entanglement generation scheme 
analyzed in RefP4 There it was shown, that all the ma- 
trix product states (MPS) of bond dimension d can be 
prepared if a string of qubits interacts sequentially with a 
(/-dimensional ancilla system and arbitrary unitaries can 
be performed on ancilla and qubit in every step. 

To apply this result to the present case, ancilla elec- 
tron and nuclear spin system together represent the con- 
trol qubit: an effective d = 2 system with Hilbert space 
spanned for given (/, m) by {\I, m, |), \I, m + 1,4)}- To 
see that arbitrary unitaries are possible, note that x- 
rotations of the control qubit are caused by resonant 
interaction, while a static Z? z -field causes z-rotations. 
From these, all single qubit gates on the control qubit 
can be constructed. The off-resonant interaction con- 
sidered before performs essentially a CNOT gate be- 
tween the passing and the control qubit. In the CNOT 
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gate the "control qubit" is the control and the passing 
electron the target, in the |T-I) an d |±) basis, respec- 
tively. Combined with single-qubit gates (on the pass- 
ing electron), this seems to be enough to allow for arbi- 
trary transformations on the coupled control-target sys- 
tem. However, the situation is more complicated since 
the effective gate performed by the off-resonant interac- 
tion differs for even and odd parity of the control qubit, 
namely V(tt) = e^^^'a^CNOT^a^, i-e., there 
is not only, as seen before, a parity-dependent phase but 
whether logical-0 or logical-1 controls the bit-flip in the 
passing qubit also depends on the parity of m. One way 
to remove this m-dependence and enable the generation 
of arbitrary states is to perform an ll I z parity measure- 
ment" by sending an electron |+) into the dot, and then 
measure it in the |±) basis after off-resonant interaction 
for a time tt. Depending on the outcome, either the odd 
or the even states are projected out. Remarkably, gaining 
this single bit of information about the 2 Ar -dimcnsional 
bath then allows us to remove all m-dependence and per- 
form clean CNOT-gates. Hence the interactions outlined 
above are sufficient to prepare all d = 2-MPS with high 
fidelity and, if the passing electrons can be brought into 
interaction with the ancilla again at any time, arbitrary 
two-qubit gates can be performed, which implies that all 
matrix product states with two dimensional bonds can 
be sequentially created^. 

Direct resonant interactions lead to very low fidelity x- 
rotations due to averaging over t he differ ent subspaces, 
indicating that prior measuremenpil22123| or coo ii n g211 f 

the spin bath might be necessary. More sophisticated 
control schemes, however, allow for near unit fidelity sin- 
gle qubit rotations with no prior preparation: In RefJ^it 
was proven that high fidelity arbitrary single qubit gates 
can be effected by a Hamiltonian H = 5<r z + Sl(<j x cos 4> + 
o y sin (j)) , when only the parameter <fi can be controlled 
precisely. For S and f2 it is sufficient to know that they 
are non-zero for some value of a controllable external pa- 
rameter and zero for another. In our situation we have 
the three Hamiltonians Hi = Aa z = B[i]a z /2 (nuclear 
Zeeman), H 2 — ■^[{m+l/2)a z +( t j l7n a x ] (resonant HF), 
and H 3 = £[±(B + m + 1/2) + B^/^ B g*)}a z (off- 
resonant HF) at hand. The Pauli matrices are acting 
on the control qubit, fii is the nuclear magnetic moment 
and £/. m = y/ljj + 1) — m(m + 1). The plus and mi- 
nus signs for H 3 can be effected through spin flips of the 
passing electron (recall that a x e lHt a x — e ia * Hcr * :t and 
<j x a z a x = —a z ). These Hamiltonians can be switched 
on and off (adiabatically 28 ) at will. Appropriate itera- 
tions of evolutions can lead to effective Hamiltonians of 
weighted sums and commutators of -Hi. 2, 3- In particu- 
lar, the subspace independence of the parameter A oc B 
allows for generation of any weighted sum of a x and a y 
with the weights being (7, m)-independent, thus making 
the results of RefP^ applicable. We have thus shown 
that while naive use of resonant interactions will lead to 
poor gate fidelities for the control qubit, enhanced control 
schemes still allow for full access to high fidelity rotations. 



Hence, in principle, universal quantum computation on 
an electron-spin quantum register can be performed, with 
all interactions mediated by the highly mixed spin bath. 

IV. EXPERIMENTAL FEASIBILITY 

We discuss now various couplings that have been ne- 
glected in the idealized Hamiltonian Eq. ([!]) but are 
present in the QD setup. We are concerned here only 
with their effects on the basic entanglement generation 
scheme. It is clear that the scheme can only work as long 
as it is fast compared to the electron T 2 time, since the 
coherence of the ancilla electron must be preserved. We 
see below that neither nuclear dynamics nor inhomogene- 
ity place more stringent conditions on our scheme. 

A. Inhomogeneity 

The HF Hamiltonian in QDs has a slightly different 
form from the one in Eq. 0, because the collective 
bath operators have a spatial dependence A/NI^ — > 
= J2i a ilii t l = i; 2 - The coupling constants at 
are oc Hi,i\ip e ( r i)\ > with |V> e (fi)| 2 being the probability 
of finding the electron at location and A — ctj 
denotes now the effective (average) hyperfine coupling 
strength. We focus our analysis on short resonant inter- 
action times At 1 <C ^/~N /A. The electronic state after 
the above protocol conditioned on a j measurement is 
proportional to 

JT( 1 /;'|A + e lK+CT|) ^ t2 +e lK+CT|)j4it M + |V')|++)xH.c., 

where the Pauli matrices act on the off-resonant elec- 
trons and \ip) — \ii . . . i^) , \ip') — \i[ . . . i' N ) label the or- 
thonormal basis of 7J eigenstates. Evaluating the matrix 
elements and introducing the normalization we get 

Pih) = E a l E (I +0 +0) + I ~j "i» x H - c -' 

1 2> 3 h,...,i N =±l/2 
if=-l/2 

/ (8) 

with the states | + (t 2 )) = e l "° t2/2 |T) + e^ * 2 / 2 ]!) and 
\~j(t 2 )) = e^^lD+e-^'* 2 / 2 !!), both of which depend 
on the nuclear spin configuration {i} via the the frequen- 
cies wo = u>o({i}) = Y^,e a iH an d = + otj and the 
normalization Affo) — J2j + cosfafa)}- The time 
dependence of the states has been omitted for brevity 
in the above formula. Straightforwardly, one now deter- 
mines the fidelity Ffa) = (4>-\p(t2)\<f>-) with the desired 
maximally entangled state ~ |+o+o) + | — — 0) as 

F(t 2 ) = 2 a ) j a p + cos^)). (9) 

3 

This expression readily gives the fidelity for arbitrary par- 
ticle numbers and arbitrary distributions of coupling con- 
stants. For N ^> 1, the obtained value is independent of 
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particle number, and for the relevant situation of Gaus- 
sian coupling F = 0.90, 0.83, 0.78 for ID, 2D and 3D, re- 
spectively. Including the difference in magnetic moments 
for Ga and As [ 75 As: ^/,As = 1-44 , 69 Ga: ^/.Ga.i = 2.02 
(60%) and 71 Ga: Mz,Ga,2 = 2.56 (40%^j these values be- 
come F — 0.83, 0.78, 0.74, indicating that our scheme is 
not compromised by realistic inhomogeneities. 

For small inhomogeneity we find the optimal time t 2 opt ^ 
by setting the time-derivative of F zero and expanding 
the equation in terms of the deviations Cj = otj — a*, 

where a* = ir/t 2 ° pt \ Going to second order in the small 
parameters tj / a* , the ensuing quadratic equation yields 

a* = m[5 - y/1 + 24(1 - W^)]/4, with Vx = ± ^ a]. 
Plugging t^ pt ^ back into Eq. (Job and keeping terms up 
to second order we find F{ir/or) = 1 — J2j( £ j/ a *) 2 - 

B. Nuclear Zeeman Energies 

For the times considered, nuclear Zeeman energies lead 
to an important relative phase -B z /i/ji2 for each two 
terms in the sum of the conditional state given in Eq. Q . 

Considering one homogeneously coupled species of nu- 
clear spins, the state of Eq. ([6| will have an additional 
m-dependent phase. In each invariant subspace this pro- 
duces an overall phase oc B z ^jt 2 m, and a relative phase 
oc B z \iiti between the two parts of the superposition. 
This might not seem harmful, but due to the parity ef- 
fect, the sign of the phase depends on the parity of to. 
Since this phase is of order tt it could spoil the protocol. 
However, by simply waiting for an appropriate time t p 
after each of the n electrons has passed the total relevant 
phase is (—l) m nB z /j,](t2 + t p ) and with t p + t2 an integer 
multiple of 7r/(/i/B 2 ) it is again TO-independent. By the 
same procedure, the nuclear Zeeman related phase can 
be removed for single-species inhomogeneous systems. 

For systems with strongly varying nuclear magnetic 
moments fii j, the relative phase depends on "which nu- 
clear spin has flipped" and the waiting time no needs to 
be chosen such that all the relative phases are close to 
2kir, otherwise the the final fidelity may be strongly de- 
graded. For the three species in GaAs this is the case, 
e.g., for B Z (J,2 +t p ) ~ 7ir, and assuming a flat wave func- 
tion still allows for a fidelity > 0.9 with moderate over- 
head in time. In principle, one can completely cancel the 
undesired phase by removing the electrons from the QD 
and reversing the magnetic field for t p = t 2 . 

C. Bath dynamics 

The major internal dynamics of nuclear spins in QDs 
stems from the indirect hyperfin e m ediated interaction, 
and the direct dipolar interactional Both mechanisms 
lead to bi-local errors that contain spin flip terms oc 
(r^' +Y'f)IZ 17 and phase changing zz-terms oc T¥l%I*. 
The transition rates for direct and indirect interactions 



are Tdd/r^ and aia>j/Q e , respectively, where fl e is the 
electron Zeeman energy and ru = |r/- — r;|. 

The dephasing interactions oc ZfJ| lead to a relative 
phase between the terms in Eq. Q , similar to the nuclear 
Zeeman energies. The energy difference, i.e. (in a mean 
field treatment) the Zeeman splitting of a single nuclear 
spin in the field of its neighbors, is a few times r^ 4 . 
Thus we need NT dd /(r%A) < 1; given T dd /rl w 0.1ms 7 
for nearest neighbors, and A w 40 ps, this condition is 
readily fulfilled even for large dots. 

We have seen above that for each term in the mixed 
state, the qubits rotate in the equatorial plane of the 
Bloch sphere with frequencies cxj when the jth nuclear 
spin has been flipped. If this particular spin is involved 
in a spin flip due to bath dynamics, the resulting rotation 
with "wrong" frequency spoils the entanglement. The er- 
rors in the rotation angle for the term containing the flip 

of the jth spin are e 3 d i = \/T^%t) 2 i a j - a k) 2 {tT) 2 , 

and the final overall errors J^. a|e^ J JfV -a 2 . We eval- 
uate the above sums in the continuum limit for Gaus- 
sian couplings and get for the indirect flips a total er- 
ror SiTr 2 A-f 2 /(-f 2 £l), where 5{ is determined by the inte- 
grals over the coupling constants. Taking A/ CI < 1 for 
the large (> 1 T) fields that we require, we find errors 
2.4%, 2.0%, and 1.5% for ID, 2D, and 3D, respectively, 
for N — 10 4 (we define N here as the number of nuclei 
within the l/(2e)-width of the Gaussian). For the di- 
rect nuclear dipole-dipole transitions the error is of size 
^d^ 2 J3^dd/(Ar^l), where numerical evaluation of the 
"dipolar integrals" 5 d yields 0.01%, 0.8%, 5% for the same 
situation as above. This overall error is thus on the order 
of a few percent for realistic situations. 



D. Storage 

We implicitly assumed the possibility of storing the 
electrons protected from any bath. In QD structures 
this could be achieved by shuttling the electrons to a 
nuclear spin-free region or employing dynamical decou- 
pling schemes, see for example 6 . The required storage 
times of a few tens of /is should be readily achieved. 



E. Imperfect Electronic Operations 

A finite probability that an up-clcctron is wrongly de- 
tected as a down electron (or vice versa), degrades the 
final entanglement. However, as only one electron needs 
to be measured, the effect is no worse for the n-partite 
GHZ-state than for the Bell state The same goes 

for variations in the resonant interaction time t\. In con- 
trast, errors in the electron preparation and variations in 
the off-resonant interaction time t<2, since they affect each 
of the n electrons lead to a fidelity reduction that scales 
exponentially with n. Variations in t 2 must be such that 
B5t2 <C 1 with B ^> A/y/N, which makes this the most 
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stringent, but still realisticP^H, requirement for electron 
timing. 

V. SUMMARY AND CONCLUSIONS 

We have considered the Heisenberg interaction of elec- 
tron spin qubits with a long-lived spin bath in a situation 
where nothing is known about the state of the bath and 
shown that nevertheless high fidelity multipartite entan- 
glement can be created via this bath. 

The qubits do neither interact directly with each other 
nor simultaneously with the bath at any time. Our pro- 
tocol thus demonstrates that even the interaction even 
with infinite temperature systems can mediate highly co- 
herent operations and thus represent a valuable resource 



for quantum information processing that merits further 
inverstigation. 

In fact, when only one bit of information is extracted 
from the spin bath, arbitrary gates between the bath and 
the qubits are possible, and all matrix product states 
with 2D bonds can be created by sequential interaction. 

The explicit protocols we presented can be realized in 
quantum dot setups and would (in typical GaAs dots) 
allow for the creation of entanglement between two elec- 
trons on a timescale of a few [is. 
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